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ABSTRACT 

We have measured the angular and spatial clustering of 671 K < 18.40, 
R — K > 5 Extremely Red Objects (EROs) from a 0.98 deg 2 sub-region of the 
NOAO Deep Wide-Field Survey (NDWFS). Our study covers nearly 5 times the 
area and has twice the sample size of any previous ERO clustering study. The 
wide field of view and B W RIK passbands of the NDWFS allow us to place 
improved constraints on the clustering of z ~ 1 EROs. We find the angular 
clustering of EROs is slightly weaker than in previous measurements, and ^(1') = 
0.25 ± 0.05 for K < 18.40 EROs. We find no significant correlation of ERO 
spatial clustering with redshift, apparent color or absolute magnitude, although 
given the uncertainties, such correlations remain plausible. We find the spatial 
clustering of K < 18.40, R — K > 5 EROs is well approximated by a power-law, 
with r = 9.7 ± 1.1 h^Mpc in comoving coordinates. This is comparable to 
the clustering of ~ 4L* early-type galaxies at z < 1, and is consistent with the 
brightest EROs being the progenitors of the most massive ellipticals. There is 
evidence of the angular clustering of EROs decreasing with increasing apparent 
magnitude, when NDWFS measurements of ERO clustering are combined with 
those from the literature. Unless the redshift distribution of K > 20 EROs is 
very broad, the spatial clustering of EROs decreases from r$ = 9.7 ± 1.1 h^ 1 Mpc 
for K < 18.40 to r ~ 7.5 h^Mpc for K > 20 EROs. 

Subject headings: cosmology: observations — large-scale structure of universe - 
galaxies: elliptical and lenticular, cD 
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1. Introduction 

The evolution of galaxy clustering is a prediction of hierarchical models of galaxy and 
structure formation (e.g., Kauffmann et al. 1999; Benson et al. 2001; Somerville et al. 2001). 
Hierarchical models for a concordance cosmology 1 predict little or no evolution of the clus- 
tering of > L* red galaxies at z < 2. Precise measurements of galaxy clustering at z ~ 1 can 
therefore test the predictions of these models. 

Extremely Red Objects (EROs; Elston, Rieke, & Rieke 1988; McCarthy, Persson, & 
West 1992; Hu & Ridgway 1994; Dey, Spinrad, & Dickinson 1995) could be the progenitors 
of local ellipticals (e.g., Spinrad et al. 1997). Roughly 80% of K s < 18.7 EROs have spectra 
with the absorption features of old stellar populations (Yan, Thompson & Soifer 2004) and 
~ 50% of K < 22 EROs have early-type morphologies (Moriondo, Cimatti, & Daddi 2000; 
Stiavelli & Treu 2001; Moustakas et al. 2004). Some EROs contain super-massive black holes, 
as ~ 15% of EROs contain an obscured Active Galactic Nucleus which can be detected by 
deep X-ray surveys (Alexander et al. 2002; Roche, Dunlop, & Almaini 2003). A direct test 
of the relationship between z ~ 1 EROs and the most massive local ellipticals is to compare 
the spatial clustering of the two populations. 

Previous constraints on the spatial correlation function of EROs, summarized in Ta- 
ble 1, are provided by pencil-beam surveys with < 0.2 deg 2 areal coverage each. Individual 
structures comprised of EROs can have sizes comparable to the field of view of these surveys 
(e.g., Daddi et al. 2000), and small surveys do not sample representative volumes of the 
Universe for highly clustered objects (e.g., Somerville et al. 2004). At z ~ 1, the transverse 
comoving distance spanned by previous ERO studies is < 20 /i~ 1 Mpc, which is much smaller 
than the size of individual structures observed in the present-day Universe. Spatial clus- 
tering measurements derived from the angular correlation function depend on ERO redshift 
distribution models. Previous angular clustering studies were unable to verify their model 
redshift distributions, as complete spectroscopic samples of EROs were unavailable. Previ- 
ous ERO spatial clustering measurements have large uncertainties and possibly large (and 
sometimes unaccounted for) systematic errors. 

In this paper, we present a measurement of the clustering of EROs using By/RIK imag- 
ing of a 0.98 deg 2 subset of the NOAO Deep Wide-Field Survey (NDWFS). The large area 
of our study provides a more representative volume than previous studies. The B^RIK 
passbands of the NDWFS allow us to constrain the ERO redshift distribution with photo- 
metric redshifts and their uncertainties. We also use photometric redshifts to select EROs 
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as a function of luminosity and redshift. We use ERO spectroscopic redshifts to verify the 
accuracy of our photometric redshifts and we compare our estimate of the ERO redshift 
distribution with spectroscopic redshift distributions from the literature. 

The structure of the paper is as follows. In §2 we provide a brief description of the 
NDWFS imaging and catalogs from which the K < 18.40 ERO sample was selected. We 
discuss our estimates of ERO photometric redshifts, and provide a comparison of ERO 
photometric and spectroscopic redshifts in §3. The selection of the ERO sample and ERO 
number counts are discussed in §4. In §5, we describe the techniques used to measure the 
angular and spatial correlation functions. The angular and spatial clustering of EROs, as 
a function of apparent magnitude, apparent color, absolute magnitude, and redshift are 
discussed in §6. We discuss the implications of our results in §7 and summarize the paper 
in §8. 

2. The NO AO Deep Wide-Field Survey 

The NDWFS is a multiband (Bw, R, I, K) survey of two ~ 9.3 deg 2 high Galactic 
latitude fields with the CTIO 4m, KPNO 4m, and KPNO 2.1m telescopes (Jannuzi & Dey 
1999). A thorough description of the optical and K-band observing strategy and data 
reduction will be provided by Jannuzi et al. and Dey et al. (both in preparation). This 
paper utilizes 0.98 deg 2 of B W RIK data in the Bootes field. B W RI imaging and catalogs 
for the entire NDWFS Bootes field became available from the NOAO Science Archive 2 on 
22 October 2004. K-band imaging and catalogs for approximately half of the Bootes field 
are also available from this archive. 

We generated object catalogs using SExtractor 2.3.2 (Bertin & Arnouts 1996), run 
in single-image mode in a similar manner to Brown et al. (2003). At faint magnitudes, 
detections in the different bands were matched if the centroids were within 1" of each other. 
At bright magnitudes, detections in the different bands were matched if the centroids were 
within an ellipse defined using the second order moments of the light distribution of the 
object 3 . Throughout this paper we use SExtractor MAG_AUTO magnitudes (which are 
similar to Kron total magnitudes; Kron 1980), due to their small uncertainties and systematic 
errors at faint magnitudes. Our clustering measurements are not particularly sensitive to 
how we measure ERO photometry, and the clustering of EROs selected with 4" diameter 



2 http : / /www. archive .noao. edu / ndwfs / 

3 This ellipse was defined with the SExtractor parameters 2 x A_WORLD, 2 x B.WORLD, and 
THETA.WORLD. 
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aperture photometry is only marginally weaker than the clustering of EROs selected with 
MAG_AUTO photometry. 

We determined the completeness as a function of magnitude by adding artificial objects 
to copies of the data and recovering them with SExtractor. To approximate z ~ 1 galaxies, 
the artificial objects have an intrinsic profile with a full width at half maximum of 0.5", 
which was then convolved with a Moffat profile model of the seeing. The 50% completeness 
limits vary within the sample area in the ranges of 26.0 < By/ < 26.7, 24.8 < R < 25.6, 
23.6 < I < 25.2 and 18.6 < K < 18.7 4 . 

Regions surrounding saturated stars were removed from the catalog to exclude (clus- 
tered) spurious objects detected in the wings of the point spread function. We excluded 
regions where the root-mean-square of the sky noise in the i^-band data was 20% higher 
than the mean, as the depth of these regions is significantly less than the mean depth across 
the field. While it is plausible that smaller variations in the sky noise could alter the mea- 
sured clustering of the faintest EROs, our main conclusions remain unchanged if we exclude 
K > 18.15 EROs from the sample. 

We used SExtractor's star-galaxy classifier to remove objects from the galaxy catalog 
which had a stellarity of > 0.7 in 2 or more bands brighter than By/ < 23.8, R < 22.8, and 
/ < 21.4. At fainter magnitudes we do not use the star-galaxy classification and correct the 
angular correlation function for the estimated stellar contamination of the sample. We do 
not use the i^-band for star-galaxy classification as there are image quality variations across 
the i^-band image stacks. We estimated stellar contamination of the galaxy sample using 
the same technique as Brown et al. (2003), where the stellar number counts were assumed 
to be a power-law and the distribution of stellar colors does not change with magnitude at 
R > 21. The contamination of the ERO sample (§4) by stars is estimated to be ~ 2%, and 
the conclusions of this paper remain unaltered unless stellar contamination is higher than 
15%. 



3. Photometric Redshifts 

Photometric redshifts were determined for all objects with / and fT-band detections. 
We provide a brief overview of the photometric redshifts here and refer the reader to our 
earlier study of 0.3 < z < 0.9 red galaxy clustering in the NDWFS (Brown et al. 2003) for a 
more detailed description of the photometric redshift code. To model galaxy spectral energy 
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distributions (SEDs), we used PEGASE2 evolutionary synthesis models (Fioc & Rocca- 
Volmerange 1997) with exponentially declining star formation rates (r models) and z = 
ages of 12 Gyr (formation z w 4). The effect of E(B — V) — 0.04 dust reddening with 
Ry = 3.1, comparable to estimates for < z < 1 early-type galaxies (Falco et al. 1999), was 
included in the r models. In Brown et al. (2003), we used models with solar metallicity at 
z — 0, which resulted in small systematic underestimates of galaxy redshifts. Simple solar 
metallicity r models underestimate the UV luminosity of galaxies (e.g., Donas, Milliard, & 
Laget 1995), so in this work we let the metallicity of the models be a function of r. This has 
the effect of slightly increasing the UV flux of the model SEDs. We verified the accuracy 
of the photometric redshifts at z < 1 with 89 M R < — 19 — 51og/i galaxies with rest-frame 
B\y — R > 1.05 and spectroscopic redshifts. After decreasing the metallicity of the models, 
the photometric redshifts of these red galaxies did not have significant systematic errors. We 
note, however, that the UV flux in galaxies can also be increased by the presence of young 
stars or by altering the properties of the dust extinction; our approach is merely a proxy 
for correcting any systematic effects in our photometric redshifts and is not meant to be 
interpreted as justifying sub-solar metallicities in the red galaxy population. We use these 
solar and sub-solar r models throughout the remainder of the paper. Color-tracks for 2 of 
the models are shown in Figure 1. For comparison, we also show two ultra-luminous infrared 
galaxy (ULIRG) templates from Devriendt, Guiderdoni, & Sadat (1999), which have bluer 
B w — R colors than the r models at z ~ 1. 

Photometric redshifts were estimated by finding the minimum value of x 2 as a function 
of redshift, spectral type (r), and luminosity. For objects not detected in the R or B w - 
bands, we estimated the probability of a non-detection using the completeness estimates 
discussed in §2. As the model SEDs do not account for the observed width of the galaxy 
locus, we increased the photometric uncertainties for the galaxies by 0.05 magnitudes (added 
in quadrature). To improve the accuracy of the photometric redshifts, the estimated redshift 
distribution of galaxies as a function of spectral type and apparent magnitude was introduced 
as a prior. The 2dF Galaxy Redshift Survey (2dFGRS) luminosity functions for different 
spectral types (Madgwick et al. 2002), with spectral evolution given by the r- models, were 
used to estimate the redshift distributions. 

We tested the reliability of the photometric redshifts with simulated galaxies and real 
galaxies with spectroscopic redshifts. Simulated galaxies were generated using the PEGASE2 
r models. The simulated data consisted of K < 18.40 galaxies with 0.6 Gyr < r < 15 Gyr 
in the redshift range < z < 5 and luminosity range 0.01 < L* < 100. The simulated object 
photometry was scattered using the estimated uncertainties, thus mimicking what would be 
present in the real catalogs. We tested the accuracy of the photometric redshifts with a few 
spectroscopic redshifts and B^RIK photometry for EROs in the NDWFS Bootes field. 
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A comparison of our photometric and spectroscopic redshifts for EROs is shown in Fig- 
ure 2. We discuss the selection criteria for the EROs in §4. The simulated galaxies in the 
left panel of Figure 2 have la uncertainties of ~ 8%. For galaxies with SEDs similar to 
the PEGASE2 r models, our procedure should yield accurate photometric redshifts. There 
are 4 NDWFS K < 18.40 EROs with spectroscopic redshifts. As shown in the right hand 
panel of Figure 2, real EROs in the NDWFS exhibit a la scatter of ~ 20% between the 
photometric and spectroscopic redshifts. There are more outliers than would be expected if 
the r models reproduced the variety of ERO SEDs. For comparison, GOODS obtains ERO 
photometric redshifts with accuracies of ~ 5% (Mobasher et al. 2004), as they have photom- 
etry and upper limits in more bands (U'U BB4 35 VV 606 RIi 775 Zs5oJHK s ). The accuracy of 
photometric redshifts is a complex function of redshift, SED and apparent magnitude. The 
accuracy of the ERO photometric redshifts could not be extrapolated from z < 1 red galaxy 
photometric redshifts, which can have uncertainties of < 10% (e.g., Brown et al. 2003). The 
accuracy of ERO photometric redshifts can not, and should not, be extrapolated from other 
samples of galaxies, such as samples selected by apparent magnitude only or from the Hub- 
ble Deep Fields (HDFs). Though our ERO photometric redshifts can only be considered 
approximations, they provide a good estimate of the ERO redshift distribution (see §5). 



4. The Extremely Red Object Sample 

We selected EROs with the R — K > 5 criterion (e.g., Elston, Rieke, & Rieke 1988; 
Daddi et al. 2000; Roche et al. 2002), though redder color cuts are sometimes used in the 
literature (e.g., Hu & Ridgway 1994; Dey et al. 1999). We have limited the sample to 
K < 18.40 EROs, to reduce the effects of completeness variations across the survey area 
on the measured clustering. As shown in Figure 3, the percentage of EROs increases from 
~ 0% of the total galaxy counts at K < 16 to ~ 8% at K ~ 18.4. 

Contamination of the ERO sample by other galaxies could significantly alter the mea- 
sured correlation function. At the magnitude limit of our sample, the uncertainty in the 
R — K color is ~ 0.25 magnitudes. For the distribution of galaxy colors shown in Figure 3, 
and assuming Gaussian photometric uncertainties, approximately 6% of the K < 18.40 ERO 
sample is contamination by R — K < 4.75 galaxies. Even if R — K < 4.75 galaxies were com- 
pletely (and implausibly) unclustered, the amplitude of the R — K > 5 angular correlation 
function would only be decreased by 12%. Contamination by 4.75 < R — K < 5.00 galaxies 
could be as high as 22% in the K < 18.40 ERO sample. This would significantly alter our 
results if the clustering of galaxies is a very strong function of color at R — K ~ 5. However, 
as discussed in §6.2, we do not see evidence of this within our dataset. Malmquist (1920) 
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bias does increase the observed number of EROs. If we assume the ERO number counts in 
Table 2 are a good approximation of the true ERO number counts, then the contribution of 
Malmquist bias to the NDWFS counts is < 8%. This would alter the measured clustering if 
ERO angular clustering is an extremely strong function of apparent magnitude. 

We assume the bulk of our sample consists of galaxies with red stellar populations. The 
colors of dusty starbursts are predicted to differ significantly from galaxies with red stellar 
populations. As shown in Figure 4, 77% of the NDWFS ERO sample has R — I > 1.15, which 
is redder than the Devriendt, Guiderdoni, & Sadat (1999) non-evolving ULIRG templates 
shown in Figure 1. Our assumption that most K < 18.40 EROs have red stellar populations 
is also consistent with the conclusions of Yan, Thompson & Soifer (2004), who find 86% of 
K s < 18.7 EROs have the absorption features of old stellar populations. 

The final sample consists of 671 objects, of which 318 are detected in the -Bty-band and 
635 are detected in the i?-band. The K < 18.4 EROs have photometric redshifts in the 
range 0.8 < z < 3.0, with the median of the distribution at z ~ 1.18. Only 5 of the 671 
K < 18.4 EROs have photometric redshifts of z > 2. ERO number counts as a function of 
.fT-band limiting magnitude are provided in Table 2 and Figure 5, along with results from 
previous surveys. We evaluated the uncertainties of the sky surface density, for the NDWFS 
and previous work, using the method discussed by Efstathiou et al. (1991), which includes 
the contribution of large-scale structure. The contribution of clustering to the uncertainties 
is typically several times larger than uncertainties determined by Poisson statistics. For 
our K < 18.4 ERO sample, accounting for the clustering increases the la uncertainty from 
5% to 20%! We note that the uncertainties quoted by some studies do not include this 
contribution (e.g., Roche et al. 2002; Miyazaki et al. 2003; Roche, Dunlop, & Almaini 2003). 
The distribution of the ERO sample on the plane of the sky is shown in Figure 6. ERO 
surveys of ~ 0.1 deg 2 often have individual structures with sizes comparable to the field of 
view (e.g., Daddi et al. 2000). While clustering and voids are evident in Figure 6, there are 
no obvious ~ 0.5° structures or gradients in the distribution of EROs in our sample. 



5. The Correlation Function 

We determined the angular correlation function using the Landy & Szalay (1993) esti- 
mator: 

DD -2DR + RR 

U{9) = RR (1) 

where DD, DR, and RR are the number of galaxy-galaxy, galaxy-random and random- 
random pairs at angular separation 9 ± 59/2. The pair counts were determined in logarith- 
mically spaced bins between 10" and 0.7°. 
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We employed the same methodology as Brown et al. (2003) to generate random object 
catalogs, correct for the integral constraint (Groth & Peebles 1977), and estimate the covari- 
ance of the oj(6) bins (using the technique of Eisenstein & Zaldarriaga 2001). The random 
object catalog contains 100 times the number of objects as the ERO catalog, so DR and RR 
are renormalized accordingly. 

The angular correlation function was assumed to be a power-law given by 

W (fl)=a;(l / )^y 7 (2) 

where 7 is a constant. This is a good approximation of the observed galaxy spatial correlation 
function from the 2dFGRS and Sloan Digital Sky Survey (SDSS) on scales of < 10 h^Mpc 
(Norberg et al. 2001, 2002; Zehavi et al. 2002). Throughout this paper we assume 7 = 1.87, 
the approximate value of 7 for z < 0.15 red galaxies from the 2dFGRS and SDSS surveys. 
For a 7 = 1.87 power-law, the integral constraint for this study was approximately 6% of 
the amplitude of the correlation function at 1'. Pair counts and the estimate of the angular 
correlation function (including the integral constraint correction) for R — K > 5.0 and 
R — K > 5.5 EROs are presented in Table 3. 

The spatial correlation function was obtained using the Limber (1954) equation; 

A0) = L Ml fM«.*.A*)^J*/(j( ^) (3) 

where ^ is the redshift distribution without clustering, £ is the spatial correlation function 
and r(9, z, z') is the comoving distance between two objects at redshifts z and z' separated 
by angle 6 on the sky. The spatial correlation function was assumed to be a power law given 
by 

ar,z) = [r/r (z)]-\ (4) 

We estimated the redshift distribution for the sample by summing the redshift likelihood 
distributions of the individual galaxies in each subsample. Model redshift distributions for 
subsamples selected by apparent magnitude and photometric redshift are shown in Figure 7. 
While the individual photometric redshifts are not especially accurate, they do include infor- 
mation provided by the observed ERO photometry and are likely to provide a fair approxi- 
mation of the ERO redshift distribution. Redshift distribution models which only reproduce 
the apparent ERO number counts and local galaxy luminosity functions (e.g., Daddi et al. 
2001; Roche et al. 2002; Roche, Dunlop, & Almaini 2003) have fewer constraints and may 
have larger systematic errors. The estimated median redshift of the K < 18.40 EROs is 1.18, 
which is almost identical to the spectroscopic median redshift of 24 K s < 18.7 EROs from 
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Yan, Thompson & Soifer (2004). The median redshift is also similar to K s < 18.5 EROs in 
the K20 spectroscopic sample (Cimatti et al. 2002, A. Cimatti 2003, private communication). 

6. The clustering of EROs 

We measured the angular and spatial correlation functions for a series of apparent 
magnitude, apparent color, absolute magnitude, and redshift bins. A power-law of the 
form u{9) = AO 1 ' 1 was fitted to the data with 7 fixed to 1.87. Much larger imaging 
surveys, including the completed NDWFS Bootes and Cetus fields, will have sufficient area 
to accurately measure 7. When parameterizing the power-law fits, we use oo(V) instead of 
c<j(1 ) as it depends less on the assumed value of 7. Using 7 = 1.80 instead of 7 = 1.87 
increases uj(1') by ~ 10% and uj(1°) by ~ 35%. The best-fit values of tq do depend on 
the assumed value of 7, but for the NDWFS ERO sample, changing 7 from 1.87 to 1.80 
increases r by only ~ 10%. Measurements of uj(V) for EROs as a function of if-band 
limiting magnitude from our study and the literature are summarized in Table 2. Angular 
correlation functions for apparent magnitude limited samples are also plotted in Figure 8. 
Estimates of oo(Y) and r for each of the NDWFS subsamples are presented in Table 4 and 
discussed in §6.1 to §6.4. 

6.1. Clustering as a function of apparent magnitude 

The amplitude of the angular correlation function for a series of apparent magnitude 
limited samples is presented in Figure 9 and Table 2, along with estimates from the literature. 
While our K < 18.4 sample has a larger volume and more objects than previous studies, our 
uncertainties are comparable to the published uncertainties of many previous studies. This is 
due to our inclusion of the covariance when fitting a power-law to the data. Our estimates of 
the amplitude of the angular correlation function are ~ 2a lower than the smaller Ks < 18.4 
ERO samples from Daddi et al. (2000). Within our sample we do not see a significant change 
in the angular clustering amplitude with apparent magnitude, but this is not unexpected as 
we span a small range of apparent magnitudes. 

The first section of Table 4 provides an estimate of ro for EROs as a function of apparent 
limiting magnitude. While the NDWFS r values are accurate to ~ 15%, each apparent 
magnitude bin spans a large range of redshift and absolute magnitude. As r is correlated 
with luminosity in other galaxy samples (e.g., Giavalisco & Dickinson 2001; Norberg et al. 
2002; Zehavi et al. 2002; Brown et al. 2003), a correlation between r and apparent magnitude 
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might be expected. We do not observe a significant correlation within the NDWFS, but our 
uncertainties are too large to rule out such a correlation. 

Combining published ERO samples provides spatial clustering measurements over a 
broad magnitude range. However, it is not possible to directly compare the published r 
measurements of different ERO samples (Table f), as different authors use different models 
of the ERO redshift distribution. Different studies also estimate the uncertainties of the 
angular correlation function and r using different techniques. In Table 3, we present the 
NDWFS pair counts for the R - K > 5.0 and R - K > 5.5 K < 18.40 ERO angular 
correlation functions, so other researchers can apply their techniques for estimating the 
correlation function to our data. Poisson statistics underestimate the uncertainties of the 
correlation function on large-scales, where object pair counts are high and the uncertainties 
of the correlation function are dominated by large-scale structure. The uncertainties of 
clustering measurements from deep pencil-beam surveys should be larger than those of the 
NDWFS, and the large scatter of K > 19 ERO clustering measurements shown in Figure 9 
may reflect this. 

If we assume the published best-fit values of the amplitude of the angular correlation 
are correct, then the angular clustering of EROs does decrease with increasing limiting 
magnitude. We find u(V) = 0.25 ± 0.05 for K < 18.40 EROs while Roche, Dunlop, & 
Almaini (2003) find u(V) ~ 0.13 for K < 22 EROs. Unless the redshift distribution of faint 
EROs is very broad, the spatial clustering of EROs is decreasing with increasing apparent 
magnitude. Several studies to measure r ~ 10 ft, _1 Mpc for faint EROs (Roche et al. 2002; 
Daddi et al. 2003; Miyazaki et al. 2003; Roche, Dunlop, & Almaini 2003), but their model 
redshift distributions contain more high redshift objects than the GOODS K s < 20.1 ERO 
photometric redshift distribution (Moustakas et al. 2004). If cj(I') ~ 0.13 for faint EROs and 
the GOODS photometric redshifts are accurate, then r decreases from 9.7 ±1.1 h^Mpc for 
K < 18.40 to r ~ 7.5 /i^Mpc for K > 20 EROs. Red galaxies at z < 1 have a comparable 
range of r values (e.g., Norberg et al. 2002; Zehavi et al. 2002; Brown et al. 2003) and their 
spatial clustering is correlated with absolute magnitude. The current measurements of ERO 
clustering are consistent with EROs being the progenitors of local red galaxies. 



6.2. Clustering as a function of apparent color 

We present the clustering of R—K > 5.0 and R—K > 5.5 EROs in Figure 9 and Table 4. 
We find the angular and spatial clustering of R — K > 5.5 galaxies does not differ significantly 
from the remainder of the sample. Low redshift galaxies may have a bimodal distribution 
of clustering properties as a function of color (Budavari et al. 2003). This could be due to 
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the bimodal distribution of galaxy colors at low redshift, or a bimodality of the clustering 
properties of galaxies as a function of star formation rate. If the clustering is bimodal at 
all redshifts, we would not expect a correlation between clustering and color within a red 
galaxy sample. We do not see a correlation of clustering with color but a larger sample with 
improved photometric redshifts is required so accurate spatial clustering measurements can 
be performed as a function of rest-frame color. 

6.3. Clustering as a function of absolute magnitude 

The clustering of EROs as a function of absolute magnitude is presented in Table 4. 
We have determined the absolute magnitudes (without evolution corrections) of the EROs 
using the best-fit r-model SED. As shown in Figure 2, ERO photometric redshifts can 
have large uncertainties and our ERO absolute magnitudes are, at best, approximations. 
The two absolute magnitude bins are approximately volume limited samples with the same 
photometric redshift range. Both absolute magnitude bins are extremely luminous, and 
contain EROs approximately 4 times brighter than the local value of L* (M£ = —23.44 ± 
0.03; Cole et al. 2001). We do not see a significant correlation between luminosity and 
r within the sample. The correlation between galaxy luminosity and clustering is only 
seen unambiguously in z < 1 samples (e.g., 2dFGRS and SDSS) which contain a factor of 
> 10 more galaxies than the NDWFS ERO sample. A strong correlation between ERO 
luminosity and spatial clustering remains plausible, and may be detected with an analysis 
of the complete NDWFS. 

6.4. Clustering as a function of redshift 

We measured the clustering of EROs within the sample with two photometric redshift 
bins, 0.80 < z < 1.15 and 1.15 < z < 1.40. We exclude EROs beyond these redshift 
ranges, as they contribute less than 10% of the total K < 18.40 ERO number counts. The 
results are presented in Table 4 and Figure 10. We do not observe significant evolution of r 
with redshift within the ERO sample. However, our uncertainties are large and the redshift 
distributions of the two samples overlap, so they are not entirely independent. 



- 12 - 



7. Discussion 

The clustering of 0.80 < z < 1.40, K < 18.40 EROs is well approximated by a power law 
with ro = 9.6 ± 1.0 h~ 1 Mpc and 7 fixed at 1.87. As EROs are thought to be the progenitors 
of local ellipticals, it is useful to compare the clustering measurements of these populations. 
In the 2dFGRS, the spatial correlation function increases from r = 6.10 ± 0.72 /i^Mpc for 
~ 2L* red galaxies to r = 9.74 ± 1.16 h^Mpc for ~ AL* red galaxies (Norberg et al. 2002). 
Other z < 1 surveys, including the NDWFS, measure comparable spatial clustering for red 
galaxies (e.g., Willmer, da Costa, & Pellegrini 1998; Budavari et al. 2003; Brown et al. 2003). 
It is not unreasonable to assume the brightest EROs are the progenitors of the ~ AL* red 
galaxies in the local Universe, as the comoving spatial clustering of the two populations is 
comparable. However, this assumes models predicting little or no evolution of the > L* 
galaxy correlation function (e.g., Kauffmann et al. 1999; Benson et al. 2001) are valid. 

If K < 18.40 EROs are the progenitors of the most luminous local red galaxies, fainter 
EROs could the progenitors of ~ L* red galaxies. Several previous studies of fainter EROs 
find they are very strongly clustered, with r ~ 10 h~ 1 Mpc (Roche et al. 2002; Daddi et 
al. 2003; Miyazaki et al. 2003; Roche, Dunlop, & Almaini 2003). This is much stronger 
than the clustering of local L* red galaxies, where r ~ 6/i _1 Mpc (Norberg et al. 2002; 
Zehavi et al. 2002). However, the ERO spatial clustering measurements could be subject 
to large, and possibly systematic, errors. Several of these measurements use model redshift 
distributions which are primarily constrained by the local galaxy luminosity function and 
faint galaxy number counts (Daddi et al. 2001; Roche et al. 2002; Roche, Dunlop, & Almaini 
2003). Other model redshift distributions use unverified photometric redshifts (Miyazaki et 
al. 2003) or redshifts which could only be verified with galaxies other than EROs (Firth et al. 
2002). As shown in Figure 9, the angular clustering of EROs does decrease with increasing 
apparent magnitude. Unless the redshift distribution of K > 20 EROs is very broad, the 
spatial clustering of K > 20 EROs is weaker than the spatial clustering of K < 18.40 EROs. 

While K < 18.40 EROs may be the progenitors of most luminous local ellipticals, our 
spatial clustering measurement should be treated with some caution. Our ERO sample 
spans broad ranges of redshift and absolute magnitude (—27.9 < M K — 51og/i < —23.9). 
The luminosity and density evolution of EROs and red galaxies at z < 1 has not been 
accurately determined. The PEGASE r models predict ~ 0.8 magnitudes of luminosity 
evolution at z < 1, so EROs would be the progenitors to ~ 2L* red galaxies. If this were the 
case, the spatial correlation function would be decreasing with decreasing redshift, which is 
unphysical. The uncertain luminosity and density evolution of EROs limits the use of > L* 
EROs to measure the evolution of the galaxy spatial correlation function. 

Current ERO spatial clustering measurements, including our study, have large uncer- 
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tainties and may be subject to systematic errors. We will significantly reduce the random 
uncertainties of ERO clustering measurements when we analyze the entire NDWFS Bootes 
field. The uncertainties and systematic errors of our photometric redshifts will be accurately 
determined as we obtain more spectroscopic redshifts. We will also improve the photomet- 
ric redshifts for EROs by using NDWFS, FLAMINGOS (Gonzalez et al. 2004), and Spitzer 
Space Telescope data. We will then be able to accurately measure the spatial clustering of 
EROs as a function of luminosity, color, and redshift. 

8. Summary 

We have measured the clustering of 671 K < 18.40 EROs with 0.98 deg 2 subset of the 
NDWFS. This study covers an area nearly 5 times larger and has twice the sample size of 
any previous ERO clustering study. The angular clustering of K < 18.40, R — K > 5.0 
EROs is well described by a power-law with cj(I') = 0.25 ±0.05 and 7 = 1.87. Using a model 
of the ERO redshift distribution derived from photometric redshifts, we find the spatial 
clustering of K < 18.40 EROs is given by r = 9.7 ± 1.0 h^Mpc comoving. Within our 
study, we detect no significant correlations between ERO clustering and apparent magnitude, 
apparent color, absolute magnitude, or redshift. However, our uncertainties are large and 
such correlations may exist. When combined with data from other studies, there is evidence 
of the angular clustering of EROs decreasing with increasing apparent magnitude. Unless 
the redshift distribution of K > 20 EROs is very broad, the spatial clustering of EROs 
decreases with increasing apparent magnitude. As the uncertainties and systematic errors 
of current ERO spatial clustering measurements are large, they do not yet provide strong 
tests of models of structure evolution and galaxy formation. 
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Fig. 1. — Color-color diagrams of the PEGASE2 r-models and non-evolving templates of 
Arp 220 and IRAS 05189-2524 (Devriendt, Guiderdoni, & Sadat 1999). The R - K = 5 
selection criterion for EROs is shown in the right-hand panel. Dots mark z = 0, 1, 2, and 3 
on the model tracks. For clarity, the redshift range shown is restricted to < z < 3 and we 
have not plotted r-models bluer than the ERO color cut. 
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Fig. 2. — A comparison of K < 18.40 R — K > 5.0 ERO photometric and spectroscopic 
redshifts for simulated and real data. Dotted diagonal lines show the measured ±lcr uncer- 
tainties of the photometric redshifts. Simulated EROs, generated using the r-models with 
photometric noise added, are shown in the left-hand panel. On the right are real EROs with 
spectroscopic redshifts. For clarity only a third of the simulated galaxies are plotted. The 
measured la uncertainties of the photometric redshifts for real EROs is ~ 20%. 
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Fig. 3. — The color-magnitude diagram and the number of K < 18.40 galaxies as a function 
of R — K color. The R — K = 5 selection criterion is marked with dashed lines. Only 25% 
of all the galaxies have been plotted in the color-magnitude diagram and R — K lower limits 
are denoted by triangles. The R magnitude limits vary across the 0.98 deg 2 sub -region, as 
the i?-band imaging consists of 4 pointings. The uncertainty of the galaxy colors increases 
from ~ 0.1 magnitudes at K — 17.00 to ~ 0.25 magnitudes at K — 18.40. R — K > 5.0 
EROs comprise ~ 8% of all K < 18.4 galaxies. 
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Fig. 4. — Color-color diagrams of the NDWFS ERO sample. Bw (left panel) and R (right 
panel) non-detections are shown with triangles and large symbols denote EROs with spec- 
troscopic redshifts. Black symbols are K — 17 or brighter while paler symbols are fainter. 
The B w — R colors of most K < 18.40 EROs are redder than the ULIRG templates plotted 
in Figure 1. A broad locus of galaxies can be seen at R — I ~ 1.5, which is coincident with 
the reddest PEGASE2 r models at 1.0 < z < 1.6. The faintest objects in the sample have 
photometric uncertainties of ~ 0.25 magnitudes, so some objects with very unusual colors 
may be photometric errors. 
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Fig. 5. — The cumulative sky surface density of R — K > 5 EROs as a function of K- 
band magnitude limit. LCIRS R — H > 4 EROs have been included using the assumption 
H — K = 1. References for each survey are listed in Table 1. We have estimated the 
uncertainties for each study using the integral constraint and assumption of square fields 
of view. The sky surface density of EROs measured by the different surveys are in good 
agreement. 
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Fig. 6. — The distribution of NDWFS EROs on the plane of the sky. EROs are shown 
with dots where the grayscale is a function of apparent magnitude (black dots are K < 17) 
while dot size is inversely proportional to the photometric redshift. Masked regions such 
as subfield boundaries, saturated stars, and non-photometric i^-band data, are shown with 
gray rectangles. Clustering is evident in the plot but there are no large gradients comparable 
to the sample area. 
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Fig. 7. — The model redshift distributions of NDWFS EROs as a function of apparent 
magnitude (left panel) and photometric redshift (right panel). The histogram of Kg < 18.7 
ERO spectroscopic redshifts from Yan, Thompson & Soifer (2004) is also shown. Unlike 
most models of the redshift distribution of K > 19 EROs (e.g., Daddi et al. 2001; Roche et al. 
2002; Roche, Dunlop, & Almaini 2003), the NDWFS K < 18.40 ERO redshift distribution 
model has few objects at z > 1.5. 
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Fig. 8. — The angular correlation function of NDWFS EROs for several apparent K-band 
magnitude limited bins. Power-law fits to the data with 7 fixed at 1.87 are shown along with 
±lcr errors (dotted lines). 
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Fig. 9. — The amplitude of the ERO angular correlation function at 1' as a function of K- 
band limiting magnitude. The error bars show the published uncertainties, which may not 
include the contribution of the covariance. The 4 arcmin 2 HDF-S measurement of K < 24 
ERO clustering (not shown) is uj(1') = 0.16 ± 0.10. Combining the NDWFS with previous 
ERO clustering studies, there is evidence for the ERO angular clustering decreasing with 
increasing magnitude. Unless the redshift distribution of K > 19 EROs is much broader than 
the redshift distribution of K < 19 EROs, the spatial clustering of faint EROs is somewhat 
weaker than the spatial clustering of K < 19 EROs. 
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Fig. 10. — The spatial correlation function of the most luminous red galaxies as a function 
of redshift. The NDWFS ERO datapoints are shown along with previous r measurements 
of z < 1 red galaxies from the 2dFGRS, SDSS and NDWFS (Norberg et al. 2002; Budavari 
et al. 2003; Brown et al. 2003). If K < 18.40 EROs are the progenitors of local ~ AL* 
early-type galaxies, the data are consistent with little or no evolution of r for z < 1.4 red 
galaxies. A ACDM model for the clustering of early-type galaxies selected by stellar mass 
(Kauffmann et al. 1999) also exhibits little evolution, but has lower clustering than most of 
the samples plotted here. 



Table I. A summary of ERO spatial clustering studies. 



Survey a 


Area 


Number of 


Magnitude 


Selection 


Additional 


Measured or 


2 distribution 15 


ro comoving c 


Assumed 




(arcmin 2 ) 


Galaxies 


Range 




selection criteria 


model z range 


model 


(h^Mpc) 


value of 7 d 


INUVVro 


oozy 


0/1 


w ^ 1 a /in 


It — I\ > 




U.o ^ 2 ^ o.U 


r flOtZj 


y. / it i.u 


l.O 1 


von 


Oz 


lo 


IS ^ 1 o o 

a < iy.z 


it — i\ > 


Uusty or orjU 


U. fyo \ 2 \ 1.41y 


Spectra 


\ Z.O 


1 Q 
l.O 


K20 


52 


15 


X < 19.2 


i?-K>5 


Old stellar SED 


0.726 < z < 1.222 


Spectra 


5.5 to 16 


1.8 


NTT-WHT 


701 


400 


K < 19.2 


i?-X>5 




0.8 < z < 2.0 


PE 


13.8 ± 1.5 


1.8 


LCIRS 


744 


337 


H < 20.0 


.R- f/ > 4 




0.7 < 2 < 1.5 


PhotZ 


11.1 ± 2.0 


1.8 


LCIRS 


407 


312 


1/ < 20.5 


.R- f/ > 4 




0.7 < 2 < 1.5 


PhotZ 


7.7 ±2.4 


1.8 


Subaru 


114 


134 


K < 20.2 


« - X 5 > 5 


Dusty SF SED 


0.0 < 2 < 4.3 


PhotZ 


12 ±2 


1.8 


Subaru 


114 


143 


X < 20.2 


R - K s > 5 


Old stellar SED 


0.0 < 2 < 4.3 


PhotZ 


11 ± 1 


1.8 


ELAIS N2 


81.5 


158 


K < 20.25 


i?-K>5 




1 <^ <3 


M-DE 


12.8 ± 1.5 


1.8 


ELAIS N2 


81.5 


158 


K < 20.25 


i?-X>5 




1 < 2 < 3 


NE 


10.3 ± 1.2 


1.8 


CDF-S 


50.4 


198 


X s < 22.0 


^775 - > 3.92 




1 < 2 < 3 


M-DE 


12.5 ± 1.2 


1.8 


HDF-S 


4 


18 


K < 24.0 


I- K > 4 






PhotZ 


16.9±|;§ 


1.8 


HDF-S 


4 


39 


X < 24.0 


/-K>3.5 






PhotZ 


6.2±f;f 


1.8 


HDF-S 


4 


23 


K < 24.0 


/ - X > 3.5 


0.8 < z < 2.0 


0.8 < 2 < 2.0 


PhotZ 


9.7 ±2.0 


1.8 



a CDF-S (Roche, Dunlop, & Almaini 2003), ELAIS N2 (Roche et al. 2002), HDF-S (Daddi ct al. 2003), NTT-WHT (Daddi et al. 2001), K20 (Daddi et al. 2002), 
LCIRS (Firth ct al. 2002), Subaru (Miyazaki ct al. 2003). 

b M-DE (merging and density evolution; Roche et al. 2002), NE (no evolution; Roche et al. 2002), PE (single burst and passive evolution; Daddi et al. 2001), PhotZ 
(photometric redshifts; Firth et al. 2002; Daddi et al. 2003, ; this work) 

c Values of ro are for a f2 m = 0.3, A = 0.7 cosmology. Uncertainties are as published, and were determined using a variety of techniques. 

d For this study, changing the value of 7 from 1.87 to 1.80 increases ro by ~ 10%. 
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Table 2. A summary of ERO angular clustering studies including number counts and sky 

surface density. 



Survey a 


Area 


Number 


of EROs 


Magnitude 


Selection 


w(l') c 


Assumed 




(arcmin 2 ) 


EROs 


per deg 2b 


Range 






7 value 




This study 


NDWFS 


3529 


256 


2.6 ±0.4 x 10 2 


K < 17.90 


R- K > 5.0 


0.36 ±0.13 


1.87 


NDWFS 


3529 


421 


4.3 ±0.6 x 10 2 


K < 18.15 


R- K > 5.0 


0.23 ±0.07 


1.87 


NDWFS 


3529 


671 


6.8 ±0.9 x 10 2 


K < 18.40 


R- K > 5.0 


0.25 ±0.05 


1.87 






Previous studes ordered by limiting 


magnitude 






NTT-WHT 


701 


58 


2.9 ±1.0 x 10 2 


K s < 18.00 


R- K s > 5.0 


0.63 ±0.26 


1.8 


NTT-WHT 


701 


106 


5.4 ± 1.8 x 10 2 


K s < 18.25 


R- K s > 5.0 


0.66 ±0.13 


1.8 


NTT-WHT 


701 


158 


8.1 ± 2.5 x 10 2 


K s < 18.40 


R- K s > 5.0 


0.58 ±0.08 


1.8 


NTT-WHT 


701 


279 


1.4 ±0.4 x 10 3 


K s < 18.80 


R- K s > 5.0 


0.37 ±0.05 


1.8 


LCIRS 


744 


337 


1.6 ±0.4 x 10 3 


H < 20.0 


R - H > 4.0 


0.33 ±0.11 


1.8 


LCIRS 


744 


201 


9.7 ±2.4 x 10 2 


H < 20.0 


/-if > 3.0 


0.36 ±0.18 


1.8 


NTT-WHT 


447.5 


281 


2.3 ±0.6 x 10 3 


K s < 19.20 


R- K s > 5.0 


0.34 ± 0.04 


1.8 


Subaru 


114 


111 


3.5 ± 1.1 x 10 3 


K s < 19.2 


R- K s > 5.0 


0.29 ±0.05 


1.8 


LCIRS 


407 


312 


2.8 ±0.5 x 10 3 


H < 20.5 


R - H > 4.0 


0.17 ±0.09 


1.8 


LCIRS 


407 


170 


1.5 ±0.3 x 10 3 


H < 20.5 


I -H > 3.0 


0.20 ±0.16 


1.8 


ELAIS N2 


81.5 


73 


3.2 ± 1.3 x 10 3 


K < 19.50 


R-K > 5.0 


0.40 ±0.16 


1.8 


ELAIS N2 


81.5 


93 


4.1 ± 1.3 x 10 3 


K < 19.75 


R- K > 5.0 


0.23 ±0.09 


1.8 


ELAIS N2 


81.5 


112 


4.9 ± 1.4 x 10 3 


K < 20.00 


R- K > 5.0 


0.20 ±0.09 


1.8 


ELAIS N2 


38.7 


63 


5.9 ± 1.5 x 10 3 


K < 20.25 


R- K > 5.0 


0.10 ±0.06 


1.8 


CDF-S 


50.4 


137 


9.8 ± 2.6 x 10 3 


K s < 21.00 


/775 - K s > 3.92 


0.14 ±0.05 


1.8 


CDF-S 


50.4 


179 


1.3 ±0.4 x 10 4 


K s < 21.50 


/?75 - K s > 3.92 


0.16 ±0.05 


1.8 


CDF-S 


50.4 


198 


1.4 ±0.4 x 10 4 


K s < 22.00 


/775 - K s > 3.92 


0.13 ±0.04 


1.8 


HDF-S 


4 


18 


1.6 ±0.8 x 10 4 


K < 24.00 


I- K > 4 


0.16 ±0.10 


1.8 



a CDF-S (Roche, Dunlop, & Almaini 2003), ELAIS N2 (Roche et al. 2002), HDF-S (Daddi et al. 2003), NTT-WHT 
(Daddi ct al. 2000), LCIRS (Firth ct al. 2002), Subaru (Miyazaki ct al. 2003). 

b The sky surface density has not corrected for the contribution of Malmquist bias, la uncertainties assume Gaussian 
errors and include the contribution of the integral constraint (using the methodology of Efstathiou ct al. 1991). 

c Uncertainties for u>(l') are as published and may not include the effect of the covariance on the uncertainty estimates. 
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Table 3. The angular correlation functions of K < 18.40 EROs from 0.98 deg 2 of the 

NDWFS. 



Color selection Angular Scales u>(0) DD DR X 10 2 RR X 10 4 



R- 


K 
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5 


.0 





.0028° 


to 





.0070° 





.631 ± 0. 


204 


88 


5515 


558532 


R — 


K 
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5 


.0 
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to 
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282 ± 0. 


100 
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33602 


3346550 


R — 


K 
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.0175° 


to 





.0440° 


0. 


.198 ± 0. 
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2274 


192715 


19207198 


R — 


K 


> 


5 


.0 





.0440° 


to 





.1106° 


0. 


.069 ± 0. 
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12084 


1115798 


107760586 


R — 


K 


> 


5 


.0 
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to 





2778° 





018 ± 0. 


037 


60806 


5811770 


555035500 


R — 


K 


> 


5 


.0 





.2778° 


to 





6977° 





018 ± 0. 


028 


244670 


23592613 


2273794800 


R — 


K 


> 


5 


.5 





.0028° 


to 





0070° 


0. 


722 ± 0. 


409 


20 


1193 


122608 


R — 


K 


> 


5 


.5 





.0070° 


to 





0175° 





064 ± 0. 


177 


80 


7502 


733626 


R — 


K 


> 


5 


.5 





.0175° 


to 





.0440° 





.240 ±0. 


087 


510 


41951 


4199066 


R — 


K 


> 


5 


.5 





.0440° 


to 





1106° 


0. 


037 ± 0. 


054 


2546 


242766 


23547906 


R — 


K 


> 


5 


.5 





1106° 


to 





2778° 


0. 


018 ± 0. 


038 


13240 


1269134 


121503332 


R — 


K 


> 


5 


.5 





.2778° 


to 





.6977° 





.020 ± 0. 


.028 


53148 


5137427 
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Table 4. The angular and spatial correlation functions of EROs from 0.98 deg 2 of the 

NDWFS. 



Selection 


Photometric 


Absolute 


Apparent 


Number 


0,(1') 


Median 


ro 


criterion 


2 range 


magnitude range 


magnitude range 


of EROs 




z 


(Ti-iMpc) 










R - K > 5.0 EROs 


selected by apparent 


magnitude 








R 


- K > 5.0 


0.80-3.00 


-27.91 < M K < -24.29 


15.92 < K < 17.90 


256 


0.36 ±0.13 


1.17 


11.0 ±2.2 


R 


- K > 5.0 


0.80-3.00 


-27.91 < M K < -23.95 


15.92 < K < 18.15 


421 


0.23 ±0.07 


1.17 


9.1 ± 1.5 


R 


- K > 5.0 


0.80-3.00 


-27.91 < M K < -23.79 


15.92 < K < 18.40 


671 


0.25 ±0.05 


1.18 


9.7 ± 1.1 


R — K > 5.5 EROs selected by apparent magnitude 


R 


- K > 5.5 


0.80-3.00 


-27.91 < M K < -24.43 


15.92 < K < 17.90 


96 


0.64 ±0.32 


1.22 


13.6 ±3.7 


R 


- K > 5.5 


0.80-3.00 


-27.91 < M K < -24.25 


15.92 < K < 18.15 


180 


0.17 ±0.13 


1.22 


7.0 ±3.1 


R 


- K > 5.5 


0.80-3.00 


-27.91 < M K < -24.11 


15.92 < K < 18.40 


314 


0.23 ±0.09 


1.24 


8.6 ± 1.9 


R — K > 5.0 EROs selected by absolute magnitude 


R 


- K > 5.0 


0.80-1.25 


-26.00 < M K < -25.00 


16.77 < K < 18.40 


108 


0.59 ±0.29 


1.16 


11.4 ±3.0 


R 


- K > 5.0 


0.80-1.25 


-25.00 < M K < -24.50 


17.35 < K < 18.40 


208 


0.28 ±0.12 


1.11 


9.0 ±2.1 








R - K > 5.0 EROs 


selected by photomd 


;ric redshift 








R 


- K > 5.0 


0.80-1.15 


-26.44 < M K < -23.79 


16.28 < K < 18.40 


318 


0.40 ±0.10 


1.04 


10.3 ± 1.4 


R 


- K > 5.0 


1.15-1.40 


-26.69 < M K < -24.46 


16.27 < K < 18.40 


292 


0.35 ±0.10 


1.28 


9.2 ± 1.4 


R 


- K > 5.0 


0.80-1.40 


-26.69 < M K < -23.79 


16.27 < K < 18.40 


610 


0.27 ±0.06 


1.15 


9.6 ± 1.0 



